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ABSTRACT 

We present deep Keck spectroscopy for 17 morphologically-selected field spheroidals in the redshift 
range 1.05 < z < 1.60 in order to investigate the continuity in physical properties between the claimed 
massive compact red galaxies ("nuggets") at z ~ 2 and well-established data for massive spheroidal 
galaxies below z ~ 1. By combining Keck-based stellar velocity dispersions with HST-h'dsed sizes, we 
find that the most massive systems (M<j yn > 10 11 M Q ) grew in size over < z < 1.6 as (i^ z )-°- 75 ± - 10 
(i.e., x2 since z = 1.5) whereas intermediate mass systems (10 11 M Q > Md yn > 10 10 M Q ) did not grow 
significantly. These trends are consistent with a picture in which more massive spheroidals formed 
at higher redshift via "wetter" mergers involving greater dissipation. To examine growth under the 
favored "dry" merger hypothesis, we also examine size growth at a fixed velocity dispersion. This 
test, uniquely possible with our dynamical data, allows us to consider the effects of "progenitor bias." 
Above our completeness limit (a > 200 km s _1 ), we find size growth consistent with that inferred for 
the mass-selected sample, thus ruling out strong progenitor bias. To maintain continuity in the growth 
of massive galaxies over the past 10 Gyr, our new results imply that size evolution over 1.3 < z < 2.3, 
a period of 1.9 Gyr, must have been even more dramatic than hitherto claimed if the red sources at 
z > 2 are truly massive and compact. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution 



1. INTRODUCTION 

The observation that many red galaxies with 
large stellar masses at z ~ 2 are 3-5 times 
more compact than equivalent e lliptic a ls in the lo- 
cal U niverse (e.s .. iDaddi et all 120051: iTruiillo et all 
20071; Ivan Dokkum et all 120081: iBuitrago et all 120081: 
Dam ianov et all 12009) has been a source of much puz- 
zlement. How can an early galaxy grow primarily in 
physical size without accreting significant stellar mass as 
required if these objects are the precursors of the most 
massive ellipticals observed today? Furthermore, studies 
of the fundamental plane and other stellar population in- 
dicators do not permit substantial recent star formation 
since z ~ 2 in massive galaxies, thus precluding growth 
by accret ion of young stars or via gas-rich ( "wet" ) merg- 
ers (e.g., iTreu et all 120051 hereafter T05). Some have 
questioned the reliability of the observations, suggest- 
ing an underestimate o f physical sizes or an overesti- 
ma te of stellar masses ( Hopkin s' et all l2009bt however, 
see iCassata et al.l 120091 for a contrasting view) . Oth- 
ers have proposed size expansion driven by self-similar 
dissipationless "dry" mergers, or mass accretion from 
minor mergers (IKhochfar k Silkll200l iNaab et alll200l 
Hopk ins" et al.ll2010l and references therein). 

To verify the compact nature of distant sources and 
to track their evolution in size and mass, it is prefer- 
able to use dynamical masses Md yn from absorption line 
spectra, which do not suffer from uncertainties associ- 
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ated with the assumed initial mass function and stel- 
lar mass estimates derived f rom broad-band photome- 
try (e.g.. iMuzzin et ail 120091) . Md yn measurements are 
available for relatively l arge samples out to z ~ 1 (T05; 
Ivan der Wei et aLll2008l hereafter vdW08), suggesting a 
small but detectable difference in average size at fixed 
mass when compared to the local universe. But beyond 
z ~ 1, there is l ittle high-quality d y namic al data for 
field spheroidals. Ivan Dokkum et all (|2009f) undertook 
an heroic observation of a single z > 2 source with a stel- 
lar mass ~2 x lO n M0 and an effective radius r e = 0.8 
kpc typical of compact galaxies at z ~ 2.3. The spectrum 
has a claimed stellar velocity dispersion of a = 510lgg 5 
km s , suggesting a remarkably dense system, van 
Dokkum et al postulate the initial dissipative collapse 
at z ~ 3 of a high mass "core" but are unable to account 
for its subsequent evolution onto the z ~ 1 scaling rela- 
tions. The quantitative effect of minor mergers on the 
physical size of a galaxy involves many variables, and it 
is unclear whether such dramatic size evolution is possi- 
ble while maintaining the tightness of the f undamental 
plane and its projections (jNipoti et a l. 2009). 

Interpretation of the observed trends at fixed Md yn is 
further complicated by the so-called "progenitor bias" 
van der Wei et alll2009l ): if galaxies grow by dry merg- 
ers, the main progenitor of a present-day massive galaxy 
did not have the same mass at 2 ~ 2. Similarly, if 
galaxies become recognizable as spheroidals only above a 
certain threshold in stellar velocity dispersion <tet that 
depends on redshift, it is clear that the addition of a 
new - and less dense - population could mimic a false 
evolutionary trend. This bias can be reduced by con- 
sidering galaxy sizes at fixed a. Foremost, a changes 
very little under a variety of growth mechanisms (e.g., 
Hopk ins" et al.lr2010T ) and it is therefore a better "label" 
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than Mdyn to track the assembly his tory. Secondly, a 
is clo sely correlated with stellar age (|van der Wei et all 
2009) and therefore offers the most direct way to track 
the evolving population. 

Given there is no clear consensus in understanding the 
continuity between the galaxy population at z < 1 and 
that at z > 2, we have embarked on a campaign to mea- 
sure a and Md yrl for a large sample of field spheroidals at 
1 < z < 1.7. This has recently become practical using 
multi-object optical spectrographs equipped with deep 
depletion red-sensitive CCDs. Our goal is to extend the 
earlier work at z < 1 (T05, vdW08) to within ~1 Gyr 
of the sample of ultracompact galaxies at z ~ 2.3. In 
this first analysis, we present new results spanning the 
redshift range 1.05 < z < 1.60. 

We adopt a ACDM cosmology with (Q m , fi„, h) = 
(0.3,0.7,0.7); all magnitudes are in the AB system. A 
Chabrier IMF is assumed where necessary. 

2. SAMPLE AND OBSERVATIONS 

Our targets were selected from archival HST/ ACS 
data in the EGS (GO 10134, PI: Davis), SSA22 (GO 
9760, PI: Abraham & GO 10403, PI: Chapman), and 
GO ODS-N (PI: Giavalis co) fields. For the EGS, we used 
the IBundv et all 0006) catalog which matches CFHT 
(BRI, ICoil eTall 120041 ugriz, CFHTLS) and Palomar 
(JK S ) photometry. Photometric redshifts are supple- 
mented by spectroscopic redshifts from the DEEP2 sur- 
vey. For SSA22, we used a photometric redshift catalog 
based on Subaru (BVRIz ) and UH 2.2m (JHK S ) imag- 
ing kindly p rovided by P ICapald (l200l. In GOODS-N, 
we used the IBundv et all (|2009D catalog which matches 
ACS and Subaru K s photometry. Gala ctic extinction 
correc tions were based on the dust maps of Schlcgcl et al. 
(1998). The parent sample for spectroscopic study in 
EGS and SSA22 was defined by I - K s > 2, I < 23.5, 
and z > 1; in GOODS-N, the photometric criteria were 
F850LP - K s > 1.5 and F850LP < 23.5. All galaxies 
satisfying these criteria were visually inspected in the 
ACS images by one of us (RSE) and those with E/S0 or 
early-disk morphology retained. 

Keck I LRIS observations were made for 14 EGS and 
SSA22 targets on 2009 June 26-28 in median seeing of 
0'.'9. The 600 mm -1 grating blazed at 1 /im was used, 
providing a velocity resolution of <7i nst = 58 km s _1 
at 9000 A. The total integration times were 40.8 ks 
and 32.4 ks in the EGS and SSA22 fields, respectively. 
On 2010 April 5-6 LRIS observations were made of 7 
GOODS-N targets with 34.8 ks of integration in 0'.'8 see- 
ing. One additional GOODS-N spectrum was secured 
with Keck II DEIMOS observations on 2010 April 11-12 
using the 831 mm -1 grating. Th e LRIS data w ere re- 
duced using the code developed bv lKelsonl (|2003l) . Spec- 
tra were extracted using optimal weighting based on 
Gaussian fits to the spatial profile. Telluric absorption 
correction and relative flux calibration were provided by 
a DA star observed at matching airmass at the end of 
each night. 

3. VELOCITY DISPERSIONS AND PHOTOMETRY 

We measured stellar velocity dispersions, a, by fit- 
ti ng broadened stellar templates using the PPXF code 
of ICappellari &; Emselleml p004). The instrumental res- 
olution was measured using unblended sky lines; their 
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Fig. 1. — Keck spectra of our 17 targets ordered by redshift (con- 
tinued in Fig. [2}. Each is smoothed with a 3 pixel (2.4 A) running 
median with sky lines omitted (black) and compared to fits to 
broadened stellar templates (red). HST images are inset with a 1" 
ruler. The order of objects matches that in Table 1. 
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Fig. 2. — Continuation of Figure [T] 
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variation with wavelength was well fit by a low-order 
polynomial. The template collection comprised 348 
stars of type F0-G 9 from the Indo-US coude library 
ijValdes et alJ l2~004) with a range of metallicities and lu- 
minosities (classes III-V). We verified that including A 
star templates does not affect our measurements. For 
each galaxy, PPXF constructed an optimal template as a 
linear combination of these stellar spectra, although our 
results do not significantly differ if the best-fitting single 
template is used. To avoid systematic errors, we masked 
pixels contaminated by OH emission. Based on tests 
with the continuum filtering, sky masking threshold, and 
stellar template choices, we assigned a systematic uncer- 
tainty to each velocity dispersion, typically 5 — 10%. We 
were able to secure a reliable dispersion for 17/22 galaxies 
(see Figs. [Hand [2]). Velocity dispersions were corrected to 
a n effective ci r cular aperture of radius R e /8 as described 
in iTreu et al.1 (|1999| ): the mean correction factor is 1.13. 

Surface photo metry was measu red in the HST images 
using GALFIT (jPeng et al.ll2002fl with a PSF determined 
from a nearby isolated star. F814W imaging was used 
in EGS and SSA22, while F850LP data were adopted in 
GOODS-N. For consistency with the local SDSS sample, 
we fit de Vaucouleurs profiles and determine circularized 
radii. We also fit Sersic profiles but found that the mean 
Seersic index n is consistent with 4 (i.e., de Vaucouleurs). 
We estimate uncertainties of ~ 10% in R e based on test- 
ing the background level, simulating the recovery of syn- 
thetic de Vaucouleur profiles placed in blank sky patches, 
and comparing with the independent measurements of 
vdW08 for the T05 subsample. 

We convert the observed ACS magnitude to the rest 
B magnitude by matching the obse rved I — K s color to 
a grid of iBruzual fc Charlotl (|2003l ) single-burst models 
of varying age and metallicity. The uncertainty in this 
fc-correction is ~ 4%. Based on the optical and NIR pho- 
tometry discussed in §2, stellar masses were estimated 
using th e Bayesian stellar po pulation analysis code devel- 
oped bv lAuger et al.l ([2009). An exponentially-decaying 
star formation history was assumed. Table 1 summarizes 
the dynamical and photometric properties for our sample 
of 17 z >1 galaxies. 

Our data allow us to compare photometrically-derived 
stellar masses to dynamical masses, which we define as 
-Mdyn = 5a 2 R e /G. Overall these are in good agreement: 
(logM dyn /Af») = 0.17 ± 0.07, consistent with the differ- 
enc e between local d y namic al and stellar masses inferred 
by ICappellari et al.l (|2006[ ) and with the dark matter 
fractions (or heavier IMFs) found by independen t meth- 
ods (e.g.. lGraves fc Faberll2010t ITreu et al.ll2010D . 

4. SIZE EVOLUTION 

Size evolution has commonly been studied by com- 
paring spheroids of the same mass at different epochs. 
Although an unlikely evolutionary path for individual 
galaxies, it is the most observationally convenient ap- 
proach, particularly at z > 1.5 where dynamical mea- 
sures are scarce. We first conduct this comparison 
in § 14.11 Following our discussion in § 1 , we then examine 
size evolution at fixed velocity dispersion in § 14.21 

4.1. Size evolution at fixed mass 

We constructed a local reference sample for com- 
parison purposes using SDSS DR7 spectroscopic data 
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Fig. 3. — Size evolution of spheroids at fixed mass. Solid 
diamonds show the mean size and its uncertainty for massive 
(log Mdyn/ Mq > 11) spheroids relative to red sequence galaxies in 
SDSS of the same dynamical mass (computed in log space). Open 
diamonds refer to intermediate-mass (10 < log M t j iyn / ~Mq < 11) 
systems. The solid line shows a (1 + z )~°-' T5 ^ z0 - 10 ^ Q ^ ne more 
massive sample. Individual galaxies in the spectroscopic samples 
are shown by colored circles. Light gray points with error bars are 
based on photometric stellar masses and show mean si zes relative 
to SD SS galaxies of the same stellar mass, based on the Sli en et al.1 
(120031) relation , for s everal samples: ISaracco et al.1 1120101 . triangle). 
Mancini ct al. (2010, pentagon), and the vdW08 (circles) compila- 
tion. 

(jStoughton et al.l I2002D . Red sequence galaxies with 
0.05 < z < 0.15 were selected using a color cut (Eq. 1 of 
lYan et aD2006l ) . Velocity disper sions were correc t ed to a 
standard aperture of i? e /8 as in Uorgensen et all (|1995f ). 
The resulting mean dynamical mass-radius relation is 
R e = 2.88 (Mdyn/10 11 M ) - 55 kpc, based on de Vau- 
couleurs radii measured in SDSS and interpolated to rest 
B band0 

Figure [3] shows the size evolution inferred by compar- 
ing the present sample and the T05 E/S0 galaxies to 
this SDSS relation. The T05 and present samples are 
well-matched in morphology and rest optical colors and 
so comprise an excellent dataset for studying evolution 
over a wide redshift interval. A simple power law fit 
(1 + z )-o-75±o.io t0 the Mdyn > 10 n m sample is indi- 
cated by the solid line. This corresponds to a 40% de- 
crease in size by z = 1, and marginally slower evolution 
than that inferred by vdW08 for M dyn > 3 x 10 10 M Q 
spheroidals ((1 + z )~°- 98±0 - 11 ; corresponding to a 50% 
decrease at Z = 1). However, the difference is partly 
explained by the correction applied by vdW08 to their 
measured sizes based on simulations. In contrast, we do 
not apply any corrections to our measured sizes given 
the lack of consensus on this matter in the litera ture (see 
also lHopkins et al.l 1200911 iMancini et al.l 120101 ) and the 
results of our own simulations. 

Our (z) = 1.3 sample probes an epoch within 2 Gyr of 
the z ~ 2.3 samples whose compact sizes have motivated 
the present work. Figure 3 illustrates that, given the size 
dispersion in the dynamical sample at a given redshift 
and the difficulties of comparing our dynamical sample 
with one whose masses are likely less precise, the sharp 
drop in size seen over this short time interval may not be 

5 The relation is consistent with that of vdW08, taking into 
account the different apertures to which we normalize velocity dis- 
persions. 
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TABLE 1 

Photometric and Spectroscopic Data 
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Note. — f denotes systems with strong Balmer absorption. * denotes the DEIMOS observation. R.A. and Dec. are in the J2000 
equinox. R. e is in kpc. Velocity dispersions (in km s _1 ) are standardized to a circular aperture of radius R e /S, and uncertainties include 
the systematic component added in quadrature. Mg is the rest-frame B band absolute magnitude. / — K s and z — K s photometry has 
been interpolated to a uniform system of Bessell I and SDSS z. Dynamical mass uncertainties reflect the total uncertainty (statistical and 
systematic) in a. Stellar masses assume a Chabrier IMF. The order of galaxies (by redshift) matches Figs. [I] and [2] 



that significant. If confirmed, however, the implied size 
evolution is quite large compared to the expected accre- 
tion from mergers over the same period, which we esti- 
mate to be 40% of a typical 10 11 M Q galaxy at z = 1.30 
Thus, if minor mergers are responsible, these results im- 
ply a very high efficiency of d log R e /d log M ~ 2.6 for 
growing spheroids, just consistent with the upper end 
of estimates determ ined from merger simulations (e.g., 
iHopkins et alllMTl . 

Also striking is the different trend seen in less mas- 
sive 10 < logMdyn/M Q < 11 galaxies (open diamonds 
in Figure [3]). Although the high-z samples are small 
and include some compact examples, we find no ev- 
idence for mean size evolution over < z < 1.6, 
i.e. oc (1 + z) 02±015 . This is consistent with a pic- 
ture in which more massive galaxies formed earlier and 
from wetter mergers with more dis sipation, creating more 
compact remnant s pheroids fe.g.. lKhochfar fc Silkll2006l : 
iTruiillo et al|[200l. and at v arianc e with the model pro- 
posed bv Ivan der Wei et al.l (|2009l ) in which lower-mass 
galaxies are the most strongly affected by progenitor bias 
and display the strongest evolution. However, we caution 
that the lower mass samples may be affected by selection 
effects, since the brighter - and therefore possibly larger 
- objects may be preferentially selected given our flux 
limits. This is not a concern for the >10 n M Q sam- 
ple, where we are complete for any reasonable mass-to- 
light ratio. A characterization of the bias requires a self- 
consistent model with Monte Carlo simulations, which is 
beyond the scope of this Letter and is left for future work 
when larger samples will be available. 

4.2. Size evolution at fixed velocity dispersion 

Comparisons at fixed mass may be affected by "progen- 
itor bias." A preferred approach, when dynamical data 
are available, is to examine galaxies of the same velocity 
dispersion. This offers the two advantages discussed in 

6 Estimated using the merger rate calculator presented in 
IHopkins efaH H2009af) . 
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Fig. 4. — Distribution of effective radii at matched velocity dis- 
persion for spheroids with a > 200 km s _1 in the present sam- 
ple (red), the T05 z > 0.9 sample (blue), and SDSS (smoothed 
grayscale with la contour). The histogram (this data) and solid 
line (SDSS) in the top panel agree by construction. The right panel 
compares the distribution of R e for the matched samples. Dashed 
and dotted lines indicate the m eans in log space. For c omparison, 
the z ~ 1.7 stacked sample of Cappellari et al. (2009) is shown 
(star), along with the si ngle galaxies of Ivan Dokk um et al.l (2009, 
hexagon, z = 2.186) and lUnodera et alj ^PloT green upper limit, 
z = 1.823). The diagonal line represents M < j yn =10 11 Mq. 

§ 1 and is illustrated in Figure |4j A cut in Mdy n includes 
only the largest galaxies at a fixed a. Therefore, if galax- 
ies below some threshold ctet are missing from the high- 
redshift samples, this could mimic an evolutionary trend 
in mass-sel ected samples. According to the preferred pre- 
scription of Ivan der Wei et all (|20 09^. ctet = 233 km s -1 
at z = 1.3. Therefore we should expect to see some effect 
for our sample, even though the presence of lower a ob- 
jects in our sample already suggests that the progenitor 
bias is not as strong. 

Figure U shows this comparison in terms of the a — R e 
plane. To avoid luminosity selection biases, we consider 
only galaxies with a > 200 km s _1 , where we are 90% 
complete, based on the SDSS distribution of luminos- 
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ity at fixed a and a conservative estimate of luminos- 
ity evolution, consistent with passive evolution of an old 
stellar population formed at zt = 3. The SDSS galaxies 
(grayscale) are weighted so as to match the a distribution 
of the 0.9 < z < 1.6 sample (top panel). The right panel 
of Figure Hthen compares the size distributions of the 
local and 0.9 < z < 1.6 samples at matching a. By fit- 
ting for size evolution at fixed a we find that sizes evolve 
as (1 + 2:)~ ' 88±0 - 19 . The good agreement with the size 
evolution inferred at fixed mass rules out strong progeni- 
tor bias. A full evolutionary model with selection effects 
is needed to include objects with lower a and quantify 
progenitor bias and size evolution more accurately. This 
is left for future work with larger samples. 

5. CONCLUSIONS 

Our Keck spectra have shown the utility of securing 
individual spectroscopic and photometric measures for a 
representative sample of z > 1 massive spheroidals. By 
probing to z ~ 1.6, we are sampling velocity dispersions, 
sizes and dynamical masses within 1.2 Gyr of the puz- 
zling population of compact red galaxies at z ~ 2.3. 

Importantly, the size evolution we infer over < z < 



1.6 at fixed dynamical mass is modest: ~ x2 for the most 
massive (log Md yn /Mo > 11) examples but much smaller 
for lower mass systems. If the compact red galaxies at 
z ~ 2 — 2.3 are their precursors, they must have grown 
dramatically in size over a very short time interval. 
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